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Abstract 

For a quantum Hamiltonian H = H{X), the observability of the energies E may be 
robust (whenever all E are real at all A) or, otherwise, conditional. Using a pseudo- 
Hermitian family of A^— state chain models H = H^^\X) we discuss some generic 
properties of conditionally observable spectra. 
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1 Introduction 



At the low energies and for the sufficiently weak interactions, quantum mechanics is 
a reliable theory. A transition to the full-fledged apparatus of relativistic quantum 
field theory is, moreover, generally believed to offer its natural extension to the higher 
energies and/or to the stronger interaction forces. In between these two extremes, 
unfortunately, there exists a huge territory of open theoretical questions. The loss 
of the internal consistency of many approximative phenomenological models may be 
encountered. Their limits of validity at certain parameters are often marked by the 
loss of the reality of the measurable quantities. One of the best known illustrations 
of such a "quantum catastrophe" is the complexification of the ground-state energy 
in the solvable model of a Dirac fermion moving in a overcritical external Coulomb 
field pp. The same effect is encountered for the Klein-Gordon boson moving in an 
external scalar field [2], etc. 

A particularly popular simulation of the breakdown of quantum stability has re- 
cently been found via the study of complex, PT— symmetric local potentials admit- 
ting the completely real spectra [3J. Such an extension of the usual phenomenology 
offers multiple advantages in the context of physics (cf., e.g., its recent review |[4j). 
At the same time, its proper implementation requires a fairly complicated mathe- 
matical apparatus ^ . This is a serious technical obstacle for any easy explanation of 
the parameter-controlled transitions between the real and complex energies. In fact, 
within the class of the PT— symmetric local potentials the reality of the spectrum 
may prove extremely fragile p]. For this reason one feels inclined to turn attention 
to the constructive study of the control of the observability of the energies (i.e., of 
the reality of all the eigenvalues of H) via some simplified models specified, say, by 
some finite-dimensional matrix Hamiltonians. 

The first results in this direction were already obtained in our two recent remarks 
[TJ [8] where we discussed some phenomenological aspects of certain "first nontrivial" , 
viz., two- and three-dimensional real-matrix models, respectively. In our present 
continuation of this effort we intend to extend our attention to the whole family of 
simplified matrix models of an arbitrary (say, even) dimension = 2J. 

The key encouragement of such a project has been found in our computer-algebra- 
based paper [9] where we revealed that a transparent mathematical structure can 
emerge not only in the above-mentioned low- dimensional models H^'^^ and H^^^ but 
also in some of their specific tridiagonal generalizations H^^^ of any dimension A^. 
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For our present purposes we pick up a subset of the latter models with the even di- 
mensions N = 2J. The reason is pragmatic - one finds just purely formal differences 
between the separate even- and odd- dimensional series of the mo dels if(^) of ref. [9]. 

Our main message will be preceded by a brief summary of the state of the art in 
section [2l A more detailed explanation of the problem will then follow in section [31 
We emphasize there that several relevant properties of our family of the Hamilto- 
nians may already be observed in its first, one-parametric real-matrix two-by-two 
member H^'^\\). We point out the intimate connection between the parity-pseudo- 
Hermiticity of the model and the mechanism which makes the "observability" (i.e., 
the reality) of the energies lost during a smooth change of the A— dependent matrix 
elements. 

The next steps of our detailed analysis of chain-models H^'^\\) are performed in 
sections m and [5] giving a detailed description of the spectra at = 4 and = 6, 
respectively. In section E] some of the quantitative conclusions of this analysis are 
found extensible to all the dimensions A^ = 2J. In particular, within the "catas- 
trophic" loss-of-the-observability scenario we stress that the fiexibility offered by the 
J independent matrix elements in if^^-^) sujfices for the simulation and enumeration 
of the eligible level-confluences igniting the imminent quantum collapse. 

Section [7] is a summary where we re-emphasize an exceptional suitability of our 
present models for deductions of some universal and generic qualitative conjectures. 

2 The loss of observability in PT— symmetric ex- 
amples 

A deep physical appeal of the "PT— symmetric quantum mechanics (PTSQM, |10j ) 
lies in the variability of its definition of the inner product in the "physical" Hilbert 
space of states [ITl [T21 [13] • The idea itself is in fact not too surprising since it 
has been discovered, forgotten and rediscovered in fleld theory [Tlj, in perturbation- 
theory mathematics [151 [16] 'w^gH ^ls in nuclear physics etc [HHTTj. After its recent 
extremely successful popularization by Bender and Boettcher [3| its use helped to 
clarify also the stability of a quantum particle in many complex quantum potentials 

mm- 

A mathematical core of the PTSQM formalism lies in its work with pseudo- 
Hermitian Hamiltonians H such that if^ = V HV^^ ^ H. The intertwiner V should 
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Figure 1: Real parts of the energies in the two-state model as functions of the 
parameter t. 

be an elementary operator which is, very often, identified with the parity [3l [T6] . 
Then, the manifestly non-Hermitian Hamiltonian H can only become selfadjoint (or, 
in the language of ref. [11], quasi- Hermitian) with respect to a nontrivial, ad hoc 
metric operator <d ^ I which depends on H = H{X) and, hence, which can also vary 
with the parameter A. 

The selected G defines the inner product in the Hilbert space of states n'^phy'^'^^^ 
so that it may become singular at some "exceptional" A [22]. This is precisely what 
paved the way to the practical physical applications of the conditional reality of the 
eigenvalues in the single-particle relativistic context of Klein-Gordon equation [23] 
and of Proca equation [21] as well as in the papers on quantum anomalies [2S], on 
supersymmetric models [26] or, beyond quantum-theory context, in cosmology [27] 
and in magnetohydrodynamics [28] . 

We believe that in virtually all of these (and many other) applications the condi- 
tional character of the reality of the eigenvalues plays the key role. For this reason 
let us now perform a certain systematic clarification of the related phenomenological 
as well as purely formal possibilities. 
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Figure 2: Imaginary parts of the energies in the two-state model as functions of the 
parameter t. 

3 The loss/emergence of observability in a matrix 
example 

The feasibihty of a controlled, conditional transition to the instability characterized 
by the complex eigenvalues may be understood as a direct consequence of the non- 
Hermitian nature of the models H = H{X) for which the parameter A either stays 
inside its "allowed" physical domain V or some of the energy eigenvalues E become 
complex. Of course, the variations of A may be expected to depend on time. For 
this reason, let us change our notation and replace A by the more explicit symbol t. 

Whenever our "time" t reaches its critical value, the related complexification 
transition may be most easily visualized in the schematic two-by-two model of ref . [7] , 




At all t < 1 this matrix is real and pseudo-Hermitian with respect to certain parity 
matrix V, 

^(2)]t^p^(2)p-l^ P= 1^ ^ j . (2) 

The specific parametrization of its matrix elements has been chosen as giving closed 
formula for the two-point spectrum, E^^ = i^/t. Obviously, these energies generated 
by the Hamiltonian ([T]) remain complex (i.e., not observable) along all the negative 
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Figure 3: The four real J = 2 energies at A = B = 1 

half-axis, t E (—oo,0) (cf. Figures [1] and [2]) . On the contrary, the matrix H'^'^^t) be- 
comes manifestly Hermitian (we could say, "conventionally physical") at t G (1, oo). 
In the middle, "unconventionally physical" interval of parameters t G (0, 1), our non- 
Hermitian Hamiltonian with real energies remains, in the terminology of the review 
paper [TT], quasi-Hermitian. 

The "unconventional" choice of t G (0, 1) leaves our matrix 

tractable as a 

valid and acceptable selfadjoint representation of an observable quantity, provided 
only that our Hilbert space of states Ti^^-* is properly re-defined (cf. the extensive 
accounts of this attitude, say, in the reviews [H |12] or in the proceedings [29]). In 
this sense the message delivered by our illustrative example ([T]) may be read as 
emphasizing that even though its form may cease to be manifestly Hermitian, the 
operator H may represent a physical observable (comment [7] may be consulted for 
all the technical details). 

Even though the symbol t represents time, let us keep both the directions of the 
t— development available. Thus, under the tacit assumption that we move to the left 
along the real t— axis in Figures ([1]) and ([2]) we speak about the confluence of the real 
energy levels followed by their subsequent complexification. This convention will be 
preferred in what follows although, alternatively, we could also change it and call the 
right-ward t— development a "big-bang-like" decomplexification of the system which 
stayed unobservable at the earlier times. 

Our specific t— parametrization of the matrix elements is privileged because it 
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Figure 4: The four real J = 2 energies ai A/2 = B = 1 

mediates the smooth t— dependence of the eigenvalues (or of their squares at least). 
Thus, even when the formal Hermiticity and/or quasi-Hermiticity of the Hamiltonian 
H^'^^ breaks down, both the Figures [1] and [2] confirm that the t = 1 division line 
between the Hermitian and non-Hermitian regime is artificial and, in any conceivable 
phenomenological context, irrelevant. 

Our two latter comments are challenging: it is not obvious what can be expected 
to happen at some higher dimensions in a suitable matrix generalization of H^'^\ A 
few answers will be offered in what follows. 



4 The first nontrivial model with four levels 



The two-parametric four-by-four chain-model 

/ -3 V^-3fj \ 

-V3-3/3 -1 2VT^ 

1 V3-3/3 

\ -V3-3/:/ 3 / 



^(4) 



(3) 



is just one of the special cases of the four-parametric pseudo-Hermitian Hamiltonians 
studied in ref. [30]. The quadruplet of the eigenenergies of matrix ([3]) is obtainable 
in closed form, 



E± = ±v^, s = s± = 3p + 2a±2J3pa + a'^ + 9f3-9a 



(4) 
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Figure 5: The four real J = 2 energies at A = 2B/3 = 1 

so that the boundary of the two-dimensional domain I) (where our pseudo-Hermitian 
matrix has the real spectrum) is composed of the two curves, viz., 

P > (^minimal = , , « G (0, 1) (5) 

and 

a > aminimal = P ~ ^ ^ /3 G (0, 1) . (6) 

Once we set 

i3 = t + Bf, a = t + Af, (7) 

we may fix the auxiliary constants A = B = 1 and stay safely inside the quasi- 
Hermiticity domain V at all the sufficiently small values of t > 0. In a way illustrated 
by Figure [3] the t— dependence of the energies remains smooth also when we cross 
the separation point t = t^*^*^) = (y^ — l)/2 ^ 0.618 between the Hermitian and 
non-Hermitian regimes with t > t^*^*^) and t < t^'^^'^\ respectively. 

Once we decide to weaken the attraction between the central energies and choose, 
say, A = 2 and B = 1, the energies split in the two well-separated pairs and they re- 
main all real whenever t^'^^) > 0.3104686356, i.e., to the right from the A— dependent 
quasi-Hermiticity boundary. The resulting t— dependence of the spectrum is dis- 
played in Figure H] where the complexifications of the two different off-diagonal ma- 
trix elements are marked by the two different vertical lines. To the right of both of 
them, our Hamiltonian H^^^ becomes standard and Hermitian, to the left of both of 
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Figure 6: The four real J = 2 energies at A = B/5 = 1 



them, our matrix if^^^ remains P— pseudo-Hermitian with respect to the usual parity 



V = 
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Between the vertical lines, an anomalous parity matrix must be chosen to define the 
pseudo-Hermiticity, 

/ -1 \ 
10 

V = 

10 
^ -1 y 

In the complementary scenario we have to weaken the attraction between the pe- 
ripheral energy levels using B > A. The choice of A = 1 accompanied by the 
weakly enhanced B = 3/2 leads to the result depicted in Figure [51 Due to the 
complexification of the most strongly attracted central pair of the levels, the quasi- 
Hermiticity is lost for t < t^^^^ = 0.2761423749. This takes place safely below the 
upper boundary t^^^^^ = 0.5485837704 of the pseudo-Hermitian regime specified by 
the standard parity operator ([8]). The latter feature is fragile. At the larger B = 5 
we get t^^^^") = 0.3582575695 which is perceivably smaller than the complexification 
bound t(Q^) = 3/5 (cf. Figure E]). 
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Figure 7: Real parts of the energies in the six-state model {A = B = C = 1). 

5 The next, more complicated chain model with 
six levels 

In the six-by-six model of ref. [9], 



^(6) 
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(1- 




^3 = 



(9) 



a 



let us set 



a = t + r + Ar 



p = t + e + Bt^ 



^ = t + e + Ct^ . 



(10) 



For A = B = C = 1 we encounter a not too interesting single pseudo-Hermiticity 
line = 0.5436890127 which is, incidentally, out of the range of Figures [7] and [H 
Thus, under our present interpretation of t as a leftwards-running time, the choice 
oiA = B = C = l would represent a fine-tuned "big crunch" process. Let us add a 
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Figure 8: Imaginary parts of the energies in the six-state model [A = B = C = 1). 

comment that after a temporary reversal of our conventional arrow of time, Figure 
[7] mimics an even more interesting scenario of a "big-bang" development. From this 
point of view there exists no observable state of our schematic system at t < and 
just a single, fully degenerate state with E = emerges at t = 0. Subsequently, teh 
evolution becomes characterized by a steady repulsion of the levels. 

After we return to our leftwards- running time convention, the peripheral-weakening 
choice of a dominant C will leave both the outermost levels far away, too weakly at- 
tracted and not sufficiently participating in the overall collapsing tendency. The 
first complexification will involve only the inner quadruplet of the energies. It may 
proceed either along the two-pair-complexification pattern indicated in Figure H] (we 
choose there A = B = C/2 = lso that the complementary pseudo-Hermiticity 
boundary moved to t^^^'^ = 1/2) or along the central-complexification pattern rep- 
resented by Figures O and [H] (necessitating an increase of B) . 

The next eligible choice of dominant A will weaken the central attraction so that 
the overall loss of the quasi-Hermiticity will be caused by the simultaneous pairwise 
mergers inside the external energy triples. The explicit decision between the two 
existing possibilities of these E ^ Q mergers will be controlled by the detailed balance 
between the size of B and C - for the weaker C the pattern will resemble Figure |H 

The last possibility corresponds to the full dominance of the constant B. This 
weakens the attraction between the central and peripheral energy pairs. A character- 
istic illustration is offered by Figure where the choice oi A = B /2 = C = lis shown 
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Figure 9: Real parts of the energies in the six-state model at A = B/2 = C = 1. 

to lead to the spectrum characterized by a central-pair-attraction dominance and by 
the related loss of quasi- Hermiticity at = 0.5157267. Our last Figure [TO] then 

shows how the other, peripheral-pair-attraction dominance modifies the previous 
result. With the choice of A = 3, B = 5 and C = 1 we get the triplet of pseudo- 
Hermiticity boundaries t^^"^ = 0.5436890127, t^^"^ = 0.4693964246 and t^^"'^ = 
0.4273046236. The model loses its quasi- Hermiticity at t < t^^^) = 0.539764657. 

6 The general chain model with N = 2J levels 

One could move on and construct various sample spectra, numerically, at a number 
of the higher even dimensions = 2J. The same family of the matrix models 
as mentioned in the previous section would still suit our purpose. The three main 
user-friendly features of all these models can be seen 

• in the fact that they represent a generalization of eq. ([9]), 

• in the nontrivial fact that the same parametrization can be used. 



• in the fairly nontrivial recommendation that 

girnax) = ^ {N - u) (cf. [9]). 




a = t + + . . . + t^-' + Gnt-^ , 



n = 1, 2, . . . , J 
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Figure 10: Real parts of the energies in the six-state model ai A/?) = B /b = C = 1. 

The numerical experiments of the two preceding sections indicate that the energy 
mergers in the spectra might admit a combinatorial classification. In such a perspec- 
tive, Figure [1] is still trivial since for the mere two available energy levels at J = 1 
there exists just their single possible merger. Still, this example enables us intro- 
duce a new convention that the two energies will be subscripted by their respective 
unperturbed integer values at gi = (giving E±i(t) at J = 1, etc). 

In this notation, the pair of the Figures H] and [5] with J = 2 documents that 
among the four available energies, there exist just the two topologically nonequivalent 
mergers. In the first case the evolution in —t connects with Ei and E_i with 
E_3, i.e., in an ordered shorthand notation, {[—3, —1], [1, 3]}. In the second case one 
connects E^ with E^^ and Ei with E^i and gets the second symbol, {[—3, 3], [— 1, 1]}. 

In such an approach we ignore the transitional multiple mergers as sampled in 
Figures [3] or [71 Still, the classification of the non-intersecting pairwise connections 
between the N = 2J energy levels requires the knowledge of the number P^^-^) of 
nonequivalent connections among the levels at every even dimension 2J. For its 
evaluation let us first order the energies in a left-right symmetric string, or lattice, 
E-{2j~i), -E_{2j-3), ■ ■ ■, E2J-3, E2J-1. This indicates that our pattern of connections 
must be left-right symmetric in this visualization. 

We already know that P^^) = 1 ^nd P(^) = 2. At J = 3 and iV = 6 it is still easy 
to check that just the following three different possibilities exist, 

{[-5, -3], [-1, 1], [3, 5]} , {[-5, 5], [-3, -1], [1, 3]} , {[-5, 5], [-3, 3], [-1, 1]} , 
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i.e., P^^^ = 3. This enumeration results not only from the required left-right sym- 
metry but also from the necessary absence of intersections between the individual 
energy-connection curves. 

Elements of the sequence of the counts of the possibilities P^^-^) will be now 
generated in recurrent manner. Firstly, it is obvious that the number p(^) of all the 
possible arrangements of the mergers at a given always incorporates a contribution 
from the subset of p(^~^) merging patterns which all contain the longest possible 
merger [— (2J— 1), 2J—1] of the two outermost energies. Let us, therefore, abbreviate 
pW _ p(^-2) — qW and, in the analysis of the remaining Q*-^"^-* options, let us 
distinguish between the even and odd J. 

In the former case with J = 2K none of the connections [±(2J — 1), m] involving 
one of the outermost energies can ever cross the center of the lattice. We have to 
fix signm = sign (±(2 J — 1)) = ±1. Thus, the absolute value of the index m, 
skipping always a pair of the levels, will run over the K—p\et of the odd integers 
2J - 3 = 4ir - 3, 2J - 7 = 4:{K - 1) - 3, 2J - 11 = 4ir - 2 - 3, . . ., 4 - 3 = 1. This 
means that we may introduce an auxiliary symbol P^^^ = 1 and evaluate, recurrently, 
as a sum of K terms at any K. This rule forms our first recurrence relation, 

p(4K) _ p{iK-2) ^ 

^ pi2K-2) . p{0) _^ p(2K~A) . p{4) _^ p{2K~6) . p{8) _^ _ _ _ _^ p(2) . piAK-8) _^ p{0) . p{4K-4) _ 

In parallel, the choice of the odd J = 2L + 1 makes the second-longest connections 
[±(2J — 1),±3] a bit shorter, leaving a two-point gap in the middle of the lattice. 
This means that the sequence of all the connections [±(2J — 1), ±n] which involve the 
outermost energies will again possess L terms. We arrive at the second recurrence 
relation, 

p(4L+2) _ p(4L) ^ 

^ p{2L-2) . p{2) _^ p{2L-4) . p{6) _^ p(2L-6) . p(10) _ ^ p(2) . p(4L-6) _^ p(0) . p(4L-2) _ 

The latter two recurrences are mutually coupled. Their numerical solution is straight- 
forward, with a sample given in Table [TJ Some of its properties are really remarkable. 
For example, empirically one finds out that the first eight (!) elements of the sequence 
p{2J) coincide with certain binomial coefficients. For the next eight elements of this 
sequence, moreover, one still finds another unexpected regularity in the differences 



P(^) = - 
2 



pirn 




p{4K+2) _ 
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Table 1: Multiplicities P(^) of the merging patterns 



K 





1 


2 


3 


4 


5 


6 


7 


8 
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piiK) 


1 


2 


6 


20 


68 


234 


808 


2798 


9700 


33656 . . . 


p(AK+2) 
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3 


10 


35 


122 


426 


1484 


5167 


17974 


62498 . . . 


RiK) 
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58 


317 


1585 


7482 . . . 


R{K) _ siK) 
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(cf. Table [I]). 

7 Conclusions 

We may summarize that our detailed quantitative analysis of schematic matrix mod- 
els revealed interesting generic qualitative features of the important phenomenolog- 
ical concept of quantum instabilities. 

Firstly we saw that the very possibility of the conditional, parameter-controlled 
emergence of the quantum collapse is closely bound to the manifestly non-Hermitian 
character of the underlying Hamiltonians. Indeed, these operators only rarely enable 
us to suppress the well known robust mathematical stability of the spectra when 
they are chosen as manifestly Hermitian. 

Of course, the removal of the Hermiticity (in the narrow sense of the invariance 
with respect to the matrix transposition and complex conjugation) does not lead to 
any conflict with the postulates of Quantum Mechanics. On the contrary, it enables 
us to make the models more flexible and more amenable to a direct control of the 
mechanism of the complexification of the eigenvalues. 

In a way related to the simplicity of our examples another key merit of them can 
be seen in the nontriviality of the related metric Q ^ I which can (and does) vary 
with the parameters. As long as the flexibility of the physics is directly encoded in 6, 
one can conclude that the access to the onset and/or breakdown of the observability 
should be mediated by the selection of the "decisive" parameters. 

We did not mention many other merits of our models (like, e.g., the particular 
advantages of their PT— symmetry, etc) because the arguments in this direction may 
be found elsewhere j4j. For compensation, let us finally note that the present outline 
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of some properties of the quasi-Hermiticity domains could be understood, in some 
sense, as the first steps towards the formulation of a certain quantum analogue of 
the Thom's theory of catastrophes [3T] . 
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Table captions 

Table m Multiplicities p(^) of the merging patterns 
Figure captions 

Figure [1]. Real parts of the energies in the two-state model 
as functions of the parameter t. 

Figure [2]. Imaginary parts of the energies in the two-state 
model as functions of the parameter t. 

Figure [3]. The four real J = 2 energies at A = B = 1 

Figure [H The four real J = 2 energies at A/2 = B = 1 

Figure [51. The four real J = 2 energies at ^4 = 25/3 = 1 

Figure [6]. The four real J = 2 energies at A = B/5 = 1 

Figure [71 Real parts of the energies in the six-state model 

{A = B = C = 1). 

Figure \Sl Imaginary parts of the energies in the six-state 
model {A = B = C = 1). 

Figure [9]. Real parts of the energies in the six-state model at 

A = B/2 = C=1. 

Figure IIOL Real parts of the energies in the six-state model 

at A/3 = B/5 = C = 1. 
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